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We present an improved QCD light-cone sum rule (LCSR) calculation of the B ^ K and Ba ^ K 
form factors by including SU(3)-symmetry breaking corrections. We use recently updated K me- 
son distribution amplitudes which incorporate the complete SU(3)-breaking structure. By applying 
the method of direct integration in the complex plane, which is presented in detail, the analyt- 
ical extraction of the imaginary parts of LCSR hard-scattering amplitudes becomes unnecessary 
and therefore the complexity of the calculation is greatly reduced. The values obtained for the 
'o6 ! relevant Bf^^) ~* K form factors are as follows: /^^^^(O) = 0.36lo:°4, /^^(O) = 0.30lo:°3, and 

O ■ /Ijf (0) = 0.38 ±0.05, /s^k(O) = 0.30 ±0.05. By comparing with the B ^ tt form factors extracted 

' recently by the same method, we find the following SU(3) violation among the B — > light form 

^ : factors: /+^(0)//+JO) = 1.38l«:iJ, ftA^)/ftM = 1-151^;^^, /Ik-(0)//I.(0) 

Ph' /I.K(0)//I.(0) = 1.17trn- 

^ '. 

XXl . PACS numbers: 13.25.Hw,12.38.Lg 



CL|' 1. INTRODUCTION 

Oh' 

1^ \ The B light meson form factors are important ingredients in the analysis of semileptonic B decays, as well as of 
. nonleptonic two-body B decays, where they serve for extraction of the Cabibbo-Kobayashi-Maskawa matrix elements. 

They have been studied by light-cone sum rule (LCSR) [H methods in several papers [I, H H i, H, 0, i, i, [13, Ell, E, [H , 
,; and most recently in p^ . In this paper we want to concentrate on the flavor SU(3)-symmetry breaking corrections 
J> ' \w B ^ K and Bg ^ K form factors, closely following the method presented in [ij. In [l^, the i? — > tt form factors 
', were analyzed, and in contrast to the previous calculations with the pole mass for mt, the MS mass Tnb{^) was 
' used. This choice more naturally follows the idea of the perturbative calculation of the hard-scattering amplitudes. 
Since the sum rule calculation of Jb and Jb^ decay constants is also available in the MS scheme [l^, we are able to 
consistently perform estimation of the fBf^^K form factors in this scheme. 
1^ ' The notion of SU(3) breaking is particularly interesting in a view of discrepancies of measured values for Bj-^,-) — s- ttK 
decay widths and CP asymmetries compared to Standard Model predictions. The K transition form factors 

enter different models for calculating these decays, and according to the recent analysis one solution of these 
discrepancies is given by assuming the large SU(3)-breaking effects, either in strong phases or in amplitudes. Our 
. ^ intention is to calculate these effects in different B(^s) K form factors by using all known SU(3)-breaking corrections 
, in the parameters and in distribution amplitudes (DAs) entering the LCSR calculation. 
j_j ■ Up to now, in [§, [13, the main SU(3)-breaking effects were included by considering SU(3)-breaking in the 
\ parameters of the leading twist DAs, such as fx / fir and hk / I^-k, and by inserting = m\ 7^ at LO. In the meantime, 
the complete SU(3)-symmetry breaking corrections in the K meson DAs are known [19]. In 19], the authors complete 
the analysis of SU(3)-breaking corrections done in pol. [21I [2^. I23L [2^ . for all twist-3 and twist-4 two- and three-particle 
DAs, by including also G-parity-breaking corrections in rus — ruq. Therefore our analysis will include complete SU(3)- 
breaking effects in both kaon DAs, as well as in the hard-scattering amplitudes at LO. At next-to- leading order (NLO) 
in the hard-scattering amplitudes, the inclusion of nig and m\ effects complicates the calculation. Because of the 
complexity of mixing between twist-2 and twist-3 DAs, we were not able to perform consistent calculations with 
rUs included in the quark propagators. Therefore, in those amplitudes we also set = m\ — 0, and consistently 
use twist-2 and twist-3 two-particle kaon DAs without mass corrections. However, we analyze the kaon mass effects 
ijp' = TO^) at NLO and include them in the error estimates. More detailed discussion about this point will be given 
in Sec. 2. 



o 
00 
o 



gorand@thphys .irb.hr| 
melic@thphys.irb.hr| 



2 




FIG. 1: Diagrams corresponding to the leading-order terms in the hard-scattering amplitudes involving the two-particle (left) 
and three-particle (right) kaon DA's shown by ovals. Solid, curly, and wavy lines represent quarks, gluons, and external currents, 
respectively. In the case of the ^ K transition, s and d quarks are interchanged, and is replaced by K" . 



Since the LO hard-scattering amphtudes are already comphcated when the twist-4 and three-particle DAs are 
included, we will use the new, numerical method to calculate the sum rules. The idea is to use the analyticity of the 
integrals, and to continue them to the complex plane. The integrals are then performed over the contour in a complex 
plane, and the imaginary part is obtained numerically. The details of the method will be given below. 



To obtain the form factors fgj^, .fsK' fsK from LCSR we consider the vacuum-to- kaon correlation function 
of a weak current and a current with the B meson quantum numbers: 



for the two different b ^ s transition currents. For definiteness, we consider the — > K (sd) flavor configuration, 
and use the isospin symmetry limit, ignoring replacement of a m quark by a d quark in the penguin current. For the 
case of fs^K^ form factors we consider the Bs K^{sd) decay. This enables us to use the same correlation function, 
with s and d quarks interchanged, but in the kaon DAs one has to take care about the fact that DAs from [l^ are 
defined for the configuration in which the momentum fraction carried by the s-quark is u (i.e., ai in the three-particle 
DAs), and u = \ — u {a2 in the three-particle DAs) is the antiquark momentum fraction. Since we want to explore 
the SU(3)-breaking corrections we will keep the kaon mass {p^ = mj^) and the rris quark mass in the DAs. The light 
quark masses will be systematically neglected, except in the ratio — rnj^Kms -\- md). 

For the large virtualities of the currents above, the correlation function is dominated by the distances = Q near 
the light cone, and factorizes to the convolution of the nonperturbative, universal part (the light-cone DA) and the 
perturbative, short-distance part, the hard-scattering amplitude, as a sum of contributions of increasing twist. In 
contrast to the pion DA, where due to the G-parity odd Gegenbauer moments vanish, the lowest twist-2 DA of a kaon 
has an expansion 



where we neglect higher moments a>2- 

We calculate here contributions up to the twist-4 in the leading order (0(a^)) and up to the twist-3 in NLO, neglecting 
the three-particle contributions at this level. Schematically, the contributions are shown in Fig.l and Fig. 2. 

By using the hadronic dispersion relation in the virtuality {p -\- q)^ of the current in the B channel, we can relate 
the correlation function ([T|) to the B ^ K matrix elements. 



2. LCSR FOR ^ K FORM FACTORS 
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FIG. 2; Diagrams contributing to the hard- scattering amplitudes at 0(as). 



Inserting hadronic states with the B- meson quantum numbers between the currents in ([T]), one isolates the B- meson 
ground-state contributions for all three invariant amplitudes F {q^ , [p + q)^) , F{q'^ , {p + q)'^), and F'^{q'^, (p + q)'^) and 
using ^ and ^ obtains 



fiKi<l') + fBKi<l') = 



F,{q\M\s^) + ^F,{q\M\s^) 



47r 



The scalar B ^ K form factor is then a combination of the vector form factor ([5]) and the form factor from 



2 2 J BK 



iBKi'f)- 



(5) 

(6) 

(7) 
®, 

(8) 



In the above, i^o(i) ^nd -fo(i) represent the LO (NLO) contributions and Jb = {Bd\mbhi^^d\Q) /m\ is the B-meson 
decay constant. As usual, the quark-hadron duality is used to approximate the heavier state contribution by in- 
troducing the effective threshold parameter s^, and the ground-state contribution of B meson is enhanced by the 
Borel-transformation in the variable {p + qY M^. Completely analogous relations are valid for Bg ^ K form 
factors, with the replacement s d in ([3]) and (|4]) and by replacing by mB^, Jb by /b^, , as well as AP by M"^ 
and Sq by Sq" in ([5] -[7]). In addition, in the derivation of the above expressions for Bg, one has to take into account 
that {Bs\bij5s\0)/m'^g^ = fBs/i^b + ms). 

The calculation will be performed in the MS scheme. The B and Bg decay constants /s^,) will be calculated in 

the MS scheme using the sum rule expressions from fis'] with 0{as, m^) accuracy. 
Each form factor can be written in a form of the dispersion relation: 



u 

F(g^Mf,),s^O-i J dse~^/<nm,F{q^,s), (9) 

ml 

where now s = {p + q)'^. The leading-order parts of the LCSR for /^^, /bk + fsK^ ^^"^ fsK form factors have the 
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following forms: 
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FIG. 3: Replacing the integration intervals by the contours in the complex planes of a = u, X, X 
in the procedure of numerical integration of LCSR amplitudes. 



^ and r2 = s/ml variables 



respectively, with X 
-3 (<^?^, d>Z^, ^•s 



ai + vas, Va = daida2da35{l — ai — 0,2 — 0^3), and with the definitions of the twist-2 {^k)^ 
twist-3 (03Xi </'3X' ^?'k), and twist-4 {4>iK, i^iK, ^iK, ^AK^ ^iK, ^ak) kaon DA's from [l^. It is easy to see that 
for — Q the above expressions resemble those given in [ij] for the B ^ n form factors. Here we have also included 
a contribution from an additional G-parity breaking twist-4 three-particle DA, 'E.^k, which was first introduced in 
[23 |. Its parameter, as well as the rest of DA parameters, are taken from where the renormalon model is used for 
describing SU(3)-symmetry breaking for twist-4 DAs. For all details about the SU(3)-symmetry breaking effects in 
the kaon DAs the reader is advised to see fioj . 

The massless (m^, nis 0) NLO contributions to the LCSR expressions for B ^ K form factors, {p + 9)^), 

etc., are the same as those given in the appendix of for B tt form factors. All features of these 0{as) 
corrections are already listed in 14], and we will not repeat them in this paper. Unfortunately, we were not able to 
perform the full NLO calculation with the mass effects included. The problem appeared by inclusion of the chirally 
noninvariant piece of the s-quark propagator, being proportional to nis, in the calculation of diagrams from Fig. 2. 
Since now p^ = mj^, the IR divergences were not present, but there appeared additional UV divergences proportional 
to nis, which have clearly exhibited the mixing among different twists. We could not achieve the cancellation of such 
singularities, since obviously some additional ingredient of mass mixing among twist-2 and twist-3 contributions was 
missing. Although interesting per se, these mixing effects are nontrivial, and there are beyond a scope of this paper. 
Hence, the repercussions of the mass effects at NLO could only be analyzed by setting nis 0. We are aware that 
keeping O(m^) effects, and neglecting the same order effect of TO<,-proportional terms is not completely justified; 
therefore, we have used the result with p^ = corrections only as an estimation for the neglected mass effects at 
NLO. 



The final LCSR expressions for B^ 
as those for the p^ — shown in [1 
present them here ^. 



K form factors, with the p 



corrections included, have a similar form 



but with a more complicated structure now. Therefore we are not going to 



3. DIRECT INTEGRATION OF THE LCSR EXPRESSIONS 



The sum rule expression for the form factors ^ requires, by definition, calculation of the imaginary part of hard- 
scattering amplitudes. Complexity of the extraction of imaginary parts of sum rule amplitudes arises already at 
the LO level, as one can notice from the expressions in Appendix A. One has to be particularly careful about the 
appearance of the surface terms there. At the NLO the results are far more complicated, as one can see in Appendix 
B of The inclusion of p^ = rri]^ effects at NLO makes the calculation even more involved. 

Therefore, we would like to present here a method which completely avoids the use of explicit imaginary parts 
of hard-scattering amplitudes, allowing one to numerically calculate amplitudes of LCSRs, analytically continuing 
integrands to the complex plane. While this method was used as a check in [l3l, here we would like to emphasis its 



^ Interested readers can obtain all expressions from the authors in Mathematica [25| form. 



6 



features and possible advantages over the traditional way of calculating sum rule amplitudes, especially when one 
performs NLO calculations. 

The main idea of the method is to deform the path of integration in order to avoid poles which are located near real 
axes. Because of the Cauchy theorem, the deformation is legitimate if the integrand is an analytic function inside the 
region bounded by the original and the new path of integration. To check the analyticity of the integrand we have 
to examine its pole structure. In the case of NLO calculations, it is also necessary to examine the position of cuts in 
logarithms and dilogarithms. Fortunately, there are just a few characteristic structures which have to be investigated. 

At LO, Eqs. (|10m2p . there are two possibilities to hit the pole, when 

ml — {q + ap)^ = 0, or q ■ p + ap^ — 0, (13) 

condition is fulfilled. In above, a = u, X, X — represents fraction of momenta between and 1, over which it has to 
be integrated. For further considerations, it is convenient to introduce the notations 

q s p m,j^ 

mi ml ml m^ 

p = (1 - a) ri + ar2 - a(l - a) r3. (14) 

By using ([14]), the conditions from (fT3|) can be written as 

2 

mg(l-p) = 0, or ^(r2-ri-r3(l-2a)) = 0. (15) 

In the case of interest r2 > ri + > 0, the second condition from above cannot be fulfilled and therefore there is 
no pole for q ■ p -\- ap^ = 0. From the first condition in psp . it follows that the integrand is approaching a pole when 
p ^ 1. That happens for the real values of r2 and a in the integration range. Note that in Eqs. (|10m2ll3)) we have 
omitted an infinitesimal imaginary quantity ie which appears in Feynman propagators. Taking it into account, the 
exact position of the pole is given by the equation 1 — p — ie = 0, which means that the poles are not located on 
the real axes of r2, but slightly below. As a consequence, one can deform the r2 path of integration into the upper 
half of the r2 complex plane to avoid passing near the poles. If poles are far away from the integration path, the 
integration is numerically completely stabile. The problem remains only when the poles coincide with the end points 
of the integration. For the integration over r2, the end points are at 1 and Sq '°'/rn^. Then p becomes 

p = a + (1 — a) ri — a(l — a) ra for r2 = 1, (16) 

s s^'"> 
p = a — ^ + (1 — a) ri — a(l - a) ra for r2 — " n > 1. (17) 
m^ ml 

In both cases p can be equal to 1 in the range of integration over a. In the case (jl6p . p is equal to 1 for a = 1, which is 
the worst possible case because this pole is located at the end point of two integrations. In p7)) . p = 1 for < a < 1, 
where, due to the specific values of ri, rs, and so, a cannot be near or 1. However, in both cases it is possible to 
move away from the poles. The complete procedure is going in this way. The first step is to shift the lower limit of 
r2 (i.e., s) integration to any point between ri -\- < 7'2 < 1- That is legitimate because all integrands are real for 
r2 < 1 and we are interested only in an imaginary part of the integrand, as can be seen from ([9]). The lower limit 
ri + ra is necessary to evade the possibility to fulfill the second condition from Eq. ([T5]) . Now we move the operation 
of taking the imaginary part in ^ outside the integral. As the third step we deform the path of the r2 integration 
into the upper half of the complex r2-plane, so that all poles are away from the integration region. For the calculation 
presented here, the new integration path is the semicircle in the complex 7'2-plane; see Fig. 3. As mentioned before, 
the pole condition still can be satisfied at the end points of integration. However, since the lower end point of the 
r2 integration is now < 1, the pole condition (p = 1) cannot be fulfilled at that end point. For the upper end point 

{f2 = Sq ''^^ 1^1) situation is as presented by Ea. p?|) . Because of the fact that this pole is in the middle of the 
range of a integration, it is possible to avoid it now by deforming the contour of a-integration into the upper half of 
the complex a-plane. Here, we again deform the integration path in the shape of the semicircle, as shown in Fig. 3. 
After that, all poles are away from the integration regions and all integrals can be performed numerically without 
facing instabilities in the integration. At the end, it remains to take the imaginary part to get the final result. 

For the NLO calculation, in addition, one has to check analytical properties of appearing logarithms and polyloga- 
rithms. It happens that for the case of the interest in one of the logarithms it is impossible to avoid crossing the cut 
when both variables r2 and a are continued to the upper half of the complex space. To avoid crossing the cut, we have 
continued only r2 to the complex plane. But now, the path of the integration for the variable a will pass near the pole 
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when r2 is approaching the endpoint Sq'°'/to^. Although the problem can be cured by the variable transformation 
and a sophisticated analytical continuation, considering the precision needed for the calculation, such a sophisticated 
method is obsolete indeed. The numerical instability shows up in the third significant digit, and therefore does not 
affect the final numerical results. 



4. UPDATED PREDICTIONS FOR THE B^,^ K FORM FACTORS 

All input parameters are listed in Appendix B. It is a compilation of the most recent determination of parameters 
entering the calculation. 

The renormalization scale is given by the expression /ijs) = y^m^^^^ — m^. Therefore, for the fg^ '^ form factors 

we use fj, = 3 GeV and for fg^j^ the renormalization scale is fig — 3.4 GeV. As usual, we will check the sensitivity of 
the results on the variation of above scales and will include it in the error estimation. 

From the general LCSR expressions for the form factors, ([5][7]), one can note that the decay constant fs (and 
correspondingly fs^ for Bg ^ K decays) enters the calculation. To reduce the dependence of the form factors on the 
input parameters, we replace fs and /s^ by two-point sum rule expressions in the MS scheme from [l^ to 0{as,m^) 
accuracy and calculate them for our preferred values of parameters. 

The usual method for deriving the working region of Borel parameters and determining effective threshold param- 
eters is used. We investigate the behavior of the perturbative expansion and smallness of the continuum contribution 
(to be less then 30% of the total contribution), and require that the derivative over the Borel parameter of the ex- 
pression for a particular decay constant, which gives the sum rule for ('t^b )i does not deviate more than 0.5 — 1% 

— 2 9 

from the experimental values for those masses. We obtain the following sets of parameters: M = 5 ± 1 GeV and 
Sq = 35.6^21 GleV^ for the i?-nieson decay constant /b calculated at /i = 3 GeV, and AI^ — 6.1 ± 1.5 GeV^ and 
Sq' = 36.6;j]}g GeV^ for Jb^ calculated at /is — 3.4 GeV. Note that the calculated central values of s^'°' follow the 
naive relation Sq" — ~ m^^ — ~ 1 GeV'^. Employing these values, the resulting decay constants are 

fs ^ 214 ± 18 MeV, fs, 250 ± 20 MeV. (18) 

In the LCSR expression for the form factors some of the uncertainties are going to cancel in the ratios, and therefore 
the error intervals of the /s and fs^ input will reduce, as one can see from the following numbers, fs = 214 ± 9 MeV 
and /b^ — 250 ±11 MeV, where the calculated error intervals come from the variation of Sq and Af^^-j only. The 

dependence of the decay constants on M^^) and s^'°' appears to be mild, as shown in Figs. [4] and [5] The SU(3) 
violation among decay constants is [T^ 

^ = 1.16 ±0.05, (19) 
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FIG. 4: Dependence of /b on (a) the Borel parameter M shown for fi — 3 GeV (solid line) and fi = 6 GeV (dashed line) and 
(b) the effective threshold parameter using the central values of all other input parameters. 
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FIG. 5: Dependence of fs^ on (a) the Borel parameter M ^ shown for /i = 3.4 GeV (solid line) and ^ = 6 GeV (dashed line) 
and (h ) the effective threshold parameter Sq" using the central values of all other input parameters. 



which nicely agrees with the values obtained from the lattice calculation and by different quark models [2^ |27|, |28| . 
The method of extraction of the Borel parameters M and Ms, and the effective thresholds Sq and Sq" for fg'^'^ 



form factors is similar to the above, and it is the same as described in We require that the subleading twist-4 
terms in the LO are small, less than 3% of the LO twist-2 term, that the NLO corrections of twist-2 and twist-3 
parts are not exceeding 30% of their LO counterparts, and that the subtracted continuum remains small, which fixes 
the allowed range of M'^^y The effective threshold parameters are again fitted so that the derivative over —l/M'^^-^ 

of the expression of the complete LCSRs for a particular form factor reproduces the physical masses m?g^ ^ with a 
high accuracy of 0(0.5% — 1%) in the stability region of the sum rules. These demands provide us the following 
central values for the sum rule parameters: = 18GeV^ = 38GeV^ = 19GeV^ and s^" = 39 GeV^ The 
dependence of the form factors on the these parameters is depicted in Figs. [B] and [71 
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FIG. 6: Dependence of on (a) the Borel parameter and (b) the effective threshold parameter Sg using the central 
values of all other input parameters. 



The complete numerical analysis yields the following predictions for the vector B ^ K and Bg ^ K form factors 
at zero momentum transfer: 



/+^(0) = 0.368 ±0.011 
/+^^(0) = 0.300 ±0.007 
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FIG. 7: Dependence of fg on (a) the Borel parameter and (b) the effective threshold parameter Sq" using the central 
values of all other input parameters. 



where the central value for /+^ is calculated at ^ = 3.0 GeV, = 18.0 GeV^, = 38 GeV^, M = 5.0 GeV^, 
and = 35.6 GeV^, and for /+ ^ at fi, = 3.4 GeV, = 19.0 GeV^, = 39 GeV^, = 6.1 GeV 



2 TT-B. 



Sq ° = 36.6 



0.10 and a^^(lGeV) = 0.25. 



GeV^. The central values for the parameters of the twist-2 kaon DA are af^(lGeV) 
The last error in and (PT|) comes from the neglected O(to^) effects at NLO. 

Finally, adding all uncertainties in quadratures, and to be on the safe side, allowing that the real mass corrections 
at NLO could reduce the final result, we obtain the following values for different B ^ K form factors: 



S6+° °5 



/^^(O) -0.38 ±0.05, 



(22) 
(23) 



and for i?s 



K form factors, 



/ii^(o) = /lK(o) = o.3ot°;l]^ 

/Ik(0) = 0.30 ±0.05. 



(24) 
(25) 



Their q dependence is shown in FigslHand FiglHl where the values in the allowed LCSR kinematical regime are shown. 
The above results for the form factors are in an overall agreement with those extracted by other methods [13, [H, [2^ . 
The predictions include also the uncertainty from the inclusion of to^ effects at NLO, which is relatively large, as 
can be deduced from (|20p and (PT|) . However, one has to be aware that this error only gives us a flavor of the size 
of neglected mass corrections at NLO, since the effects in the hard-scattering amplitude could not be included, 
and we expect that there will be a partial cancellation among and n{\^ contributions, being of similar size. At the 
leading order, the inclusion of effects (appearing only in DAs) and m|. effects reduces the results by 2.5% — 4%, 
depending of the value of g^. 

In order to be able to comment on the SU(3)-breaking effects, these values have to be compared with the results 
obtained by the same method for B ^ tt form factors [l4l|, which we quote here: 



26+° °'* 



.(0) = 0.255 ±0.035 



(26) 
(27) 



By varying parameters in a correlated way, finally we predict the following SU(3)-breaking ratios: 



1 oo+O.ll 
-'^•'-'"-O.IO ' 



1 40+0-18 



/Xk(0) 



1.15 



+0.17 
-0.09 ' 



1.17] 



0.11 



(28) 
(29) 



The complete SU(3) violation comes from SU(3)-breaking corrections in all parameters, mainly from fx/ fir, M-R'/Mtt, 
from the difference in the sum rule parameters, sb and M, as well as from the difference in the /b^ and fs ratio. 
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FIG. 8: The LCSR prediction for form factors fgi^iq^) (solid line), fsKi'l'^) (dashed line), and f^jxiff') (dash-dotted line) at 
< < 12 GeV^ and for fj, = 3 GeV, Sq = 38 GeV^, = 18GeV^, and the central values of all other input parameters. 



+ , r, 
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FIG. 9: The LCSR prediction for form factors fg k{i^) (solid line), fssKiff) (dashed line), and f's^Kil 



atO<q^ <12 GeV^ and for = 3.4 GeV, Sq 



39GeV^ M: 



(dash- dotted line) 
19 GeV^, and the central values of all other input parameters. 



Compared with the values from the second paper in [l7j . where the similar LCSR analysis was done, we find nice 
agreement with the results in 

It is also interesting to explore an overall SU(3)-breaking factor, which appears in factorization models for — > 
Ktt amplitudes [H, [SO] : 



U f 



1 01+°°^ 
-■^•^^-o.is- 



(30) 



-K 



For the /k/ f-w ratio we use (jB3l) . Although there is a SU(3) violation among form factors and in the masses, the 
predicted value for ^ shows almost exact SU(3) symmetry. On the other hand, the above ratio enters the prediction for 
Bs — > K~iT^ amplitude obtained by employing U-spin symmetry. U-spin symmetry cannot be trusted [17], as we can 
note by inspecting another U-spin relation. By neglecting penguin and annihilation contributions, under the U-spin 
symmetry assumption Afact{Bs K^) / Afact{Bd tt+tt^) ~ 1, [itI. [30I. [31|| . while our prediction amounts to 



AfactiBs ^ K+K-) _ fK ftKi^K) - m 



i-fact 



{Bd 



= 1.41 



mt 



+0.20 
-0.11 ' 



(31) 



a quite substantial U-spin violation. 
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5. SUMMARY 



In this paper we have investigated the SU(3)-symmetry breaking efTects in the B K and Bg ^ K form factors. 
The analysis has involved the SU(3)-breaking corrections both in the LO (up to twist-4 corrections), as well as in 
the NLO calculation, estimating SU(3) corrections for the twist-2 and twist-3 contributions. Although at NLO we 
were not able to consistently include 0{ms) ^ 0{mj^) effects, we have included effects in the error analysis of 
our results. We have presented a method of numerical integration of sum rule amplitudes, which greatly facilitates 
the calculation, especially the calculation of the radiative corrections. By investigating some of the SU(3) and U-spin 
symmetry relations, we have shown that such relations have to be considered with a precaution, since some of them 
can be badly broken. 
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APPENDIX A: EXPLICIT FORMULAS FOR THE LEADING-ORDER LCSR EXPRESSIONS 

Although the intention of this paper is to promote the numerical method for calculating LCSR amplitudes, for 
which the explicit expressions for the imaginary parts are superfluous, because the result can be obtained by direct 
integration of the starting amplitudes (at LO they are given by Eqs. (|10lll2|) ). we have decided to list here the LO 
LCSR expressions for fs^^^jK form factors, since to our best knowledge, these expressions were never clearly presented 
in a form which includes complete mass corrections. 

The LO part of the /^^^ LCSR, (O, has the following form: 



Fo{q',M^,s^)^mifK I due 



rrib 



m; 



1 



Ull)iK{u) + 1 



2u^m|- 



(m 



.2 2 ■sd<p3Kiu) 



2 2 \ 



du 



dvip4:K{v) 



du^ 



12um^ 



mi. — + u^mj^ du (m^ — + u^m^^Y 



d 

du 



2um^pr 



K 



du 



2ura\ 



— q^ + u^m^ 



U—+11 

du 



iu'^mj^ 



rri^ — q"^ + u^m^ 



^ Iak{.u) 



2um?j^{m1 



? 2 2 

q — u m 



K) 



[ml — q^ + u'^mj^)^ 

2 J- --^ 

my KG ^ 



d 

du 

lWW 



6um'^ 



dw^rr(«) + 



4 mi — cp 



'if 



dK,^{u) 
du 



t>iK\U) 



(Al) 



where ' 



-u, Uo = 



I)) /(2m2,), 



and the short-hand notations 
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introduced for the integrals over three-particle DA's are 



hK{u) 



hK{u) 



u i 

j dai j [Avp- q- {l-2v)um\]^zK{a^) 



(«-Qi)/(l-Qi) 
1 

dai 

u 1 



Q!2 = 1 - ai - (33, 

as = (u - ai)/t; 



V 



da I 

M 1 



dv 

V 



a2 = 1 — ai — as, 
0:3 = (u — ai)/v 



a2 = 1 — ai — a3, 
a3 = (u - 







(if 



w(l - v)^4K{ai) 



^2 = 1 - «! — as, 
as = - ai)/w 



(A2) 



The twist-4 two-particle DAs are defined with the help of the two DAs <p4K = 'PJk+^I'Yk^ and iI^4k = -^PIk+^Vk^ ■ 
The leading-order LCSR for /^^, + fgj^, looks like 



1 mfc \ u 



1 



i>4K{u) 



2urrij^ 



+ 



du 7 



^ + u^rnj^ \ du^ mt 



hK{u) 
d 



L 6u du 

dvij^iK {v) 



q^ + u^m^pr du (to^ — g'^ + v?m\) 



diUKii) 



dv^fj^{v) 



where 



(A3) 



hK{u) = da 



dv 



[(3-2i.)]$3K(a,) 



(u-ai)/(l-ai) 

Finally, the leading-order LCSR for the penguin form factor, reads 



a2 = 1 - ai - as, 
as = (u - ai)/v 



(A4) 



Fo^(<;^M^s^) = mb/K / due 



3(m^ — (7^ + u'^mj^) 
1 



du 



— 2,2 2 '^3K (w) 



li — q^ + u^m 



2um?p- 



du rri^ — q^ + u'^rnj^ 



6umj^ 



du 



12um\ 



4(to^ — + u^rn^j^) \du'^ — q^ + u^ra^ du (m^ — q^ + u'^rnj^)^ 



(t>AK{u) 



ml - q'^ + m|. 



4 m? — q'^ + m 



K 



d4>7irin) 

du 



(A5) 
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Parameter 


Value at /i = 1 GeV 




128 ± 21 MeV 


K 


0.10 ±0.04 [35] 


af 


0.25 ±0.15 







fsK 


0.0045 ± 0.0015 GeV^ 




-1.2 ±0.7 




1.6 ±0.4 


S'k 


0.20 ± 0.06 GeV^ 


UJ4K 


0.2 ±0.1 




-0.09 ±0.02 



TABLE I: Input parameters for the kaon DA 's fl2 . [71/. 



Parameter 


Value 




0.1176 ±0.002 


{qq){l GeV) 
(ss)/{qq) 

mo 


-{2A6t\l MeV)3 
0.8 ±0.3 

0.012in?2 GeV^ 
0.8±0.2GeV^ 



TABLE IL Additional input parameters for the /b and fs^ sum rules. 



and 



U i 



2*4/f(a.)-(l-2t;)$4K(aO 



(«-Qi)/(l-Ql) 

±2(l-2w)$4K(aO-$4K(az) 



a2 = 1 - ai - as, 
"3 = (w - C(l)/V 



(A6) 



Note the appearance of the surface terms in form factors above. Being proportional to the Wandzura-Wilczek part 
of (jjiK and ip4K, they vanish for — > 0. 

The expressions for f^'^j^ form factors follows from above, by replacing u by 1 — u and by interchanging ai and 

a2, and rrid and rris in the kaon DAs, i.e., by replacing DAs of by those for K^. 



APPENDIX B: PARAMETERS USED IN THE CALCULATION 



In this appendix we summarize the parameters used in the calculation of fsf^^jK form factors as well as in the 

calculation of Jb^^^ Tables HI - IIIII For DAs, the parameters and their /i dependence are taken from [U, [l^. The MS 
mass rnb entering the calculation is (s^ 

TOb(mb) =4.164 ± 0.025 GeV. (Bl) 

The as{mz) value is the Particle Data Group (PDG) average [33|. There is a new value for the fx decay constant 
[H, prepared for the PDG's 2008 edition, 

= (156 ±0.2 ±0.8 ±0.2) MeV, (B2) 
which central value we adopt here, and the ratio of and tt^ decay constants is given by 

^ = 1.196 ±0.002 ±0.006 ±0.001. (B3) 
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90 


= 


IK 
9o = 


1 x2 


^-0^ = 


1 ^2 




= 


1 r2 K 

"0 — -^"ko.1 


0f 


- T2"K 


IK 

91 = 




V'f = 


7 r2 




- 7 ^2 




92 


— ^20"Kai 


IK 
92 = 




V'2^ = 




^2^ 







TABLE III: Three-particle twist-4 parameter relations for kaon DA 's derived from the renormalon model fldl. IKj l. 



The rfis mass is the average [35| of the QCD sum rule determinations from [36|, |37| and covers the rris mass range 
given the most recently in [3^ . For the twist-2 kaon DA we use the first Gegenbauer moment calculated at NNLO 
accuracy from [sll . Since the existing fits and calculations for the value of the second Gegenbauer moment af^ show 
small SU(3) violation, with the large error, we accept here that af = aj, with the value for af^ given in Table 1 p^ . 
We also use rriK = 497.648 ± 0.022 MeV, uib = 5279.5 ± 0.5 MeV, and ms, = 5366.1 ± 0.6 MeV t33|]. 
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